Interferometric approach to measuring band topology in 2D optical lattices 
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Recently, optical lattices with non-zero Berry's phases of Bloch bands have been realized. New 
approaches for measuring Berry's phases and topological properties of bands with experimental tools 
appropriate for ultracold atoms need to be developed. In this paper, we propose an interferometric 
method for measuring Berry's phases of two dimensional Bloch bands. The key idea is to use a 
combination of Ramsey interference and Bloch oscillations to measure Zak phases, i.e. Berry's 
phases for closed trajectories corresponding to reciprocal lattice vectors. We demonstrate that this 
technique can be used to measure Berry curvature of Bloch bands, the tt Berry's phase of Dirac 
points, and the first Chern number of topological bands. We discuss several experimentally feasible 
realizations of this technique, which make it robust against low-frequency magnetic noise. 
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Introduction. Topology underlies many fundamental 
physical phenomena in two-dimensional materials, most 
notably the quantum and anomalous Hall effects [iHl]- 
The topological features of a 2D Bloch band are de- 
termined by the Berry's (geometrical) phases [5] - the 
phases picked up during an adiabatic motion of a particle 
along closed trajectories in quasi-momentum space. The 
energy bands are classified by an integer-valued topolog- 
ical invariant C, the first Chern number, which is propor- 
tional to the Berry's phase for a trajectory enclosing a 
full Brillouin zone (BZ). A filled band with Chern num- 
ber C is characterized by the quantized Hall conductitivy 
CeV/i [HE]. 

Over the past few years, the interest in topological 
properties of 2D systems was strongly revived, following 
the discovery of new materials, including graphene [6 
and topological insulators |7], where the Berry's phases 
play an important role in defining the transport proper- 
ties. For example, the band structure of graphene hosts 
two massless Dirac points at the corners of the BZ. A 
trajectory enclosing a Dirac point is characterized by the 
Berry's phase tt, and it is this tt Berry's phase that lies 
at the heart of new phenomena observed in graphene, 
such as the half-integer quantum Hall effect and weak 
anti-localization^ . 

There is currently a strong interest in realizing topo- 
logical band structures in cold atomic systems [8til3j. 
Very recently, highly tunable 2D optical lattices with 
non-trivial Berry's phases have been demonstrated exper- 
imentally [HHH]. In particular, staggered flux lattice [5], 
as well as brick-wall and honeycomb optical lattices with 
massless Dirac points have been engineered [SlITU]. More- 
over, several intriguing theoretical proposals of how to 
realize topologically non-trivial optical lattices with non- 
zero Chern number were put forward |14H23j . 

However, the traditional transport measurements, 
which can be used to determine the Chern number of 




FIG. 1: Experimental setup for measuring the Zak phase. A 
cloud of ultracold atoms with a well-defined quasi-momentum 
ko is loaded into a 2D optical lattice. Gi,G2 denote recip- 
rocal lattice vectors. Initially, atoms have spin up. (a) A 
7r/2 pulse creates a coherent superposition of j t) and | l) 
states (marked by blue and red). After that, spin-selective 
forces ±F parallel to Gi are applied, (b) After half a pe- 
riod of Bloch oscillations, when the two spins meet in the 
quasi-momentum space, another 7r/2 pulse is applied. The 
accumulated phase difference between the two states, which 
contains the Zak phase contribution, is measured by reading 
out the phase of the resulting Ramsey fringe. 



a band via the quantized Hall effect, are very challenging 
in cold atomic systems. Thus, novel methods of probing 
the topology of bands in optical lattice are needed. One 
possible route to studying the topological structure of op- 
tical lattices relies on the fact that the Berry curvature 
gives rise to an anomalous velocity in semiclassical dy- 
namics that can be monitored through in-situ images of 
the atom cloud [2] [24] . Another possibility that has been 
pointed out recently is that time-of-flight images could 
be used to reveal topological invariants [551 ES] ■ 
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Here we propose an alternative, interferometric 
method for measuring Berry's phases, Berry curvature 
and Chern number of bands in 2D optical lattices. The 
idea is to combine Bloch oscillations with Ramsey inter- 
ferometry for particles with two internal states, |t), 
loaded into the optical lattice. The BZ has the topology 
of a torus, and therefore during the Bloch oscillations 
particles follow closed trajectories corresponding to the 
cycles of the BZ torus. The Berry's phases of such tra- 
jectories are known as the Zak phases |27) . 

We note that the two main ingredients required for 
the implementation of our proposal - Bloch oscillations 
and Ramsey interferometry, have become standard tech- 
niques in cold atoms experiments [281 - 133] : thus, we ex- 
pect that our proposal can be used to measure topolog- 
ical properties of 2D lattices in the near future. Very 
recently we have successfully applied a related approach 
to measure the Zak phases of topological Bloch bands in 
one dimension |34j . 

We show how the measurements of the Zak phase allow 
one to obtain the more familiar topological characteris- 
tics of 2D Bloch bands. First, measuring the change of 
the Zak phase in the BZ, one can determine the distribu- 
tion of the Berry curvature. This allows one to measure 
the Chern number of the band, given by the winding 
number of the Zak phase across BZ. Second, we show 
that the difference of Zak phases measured along certain 
trajectories can be linked to the Berry's phase of Dirac 
fermions. To be concrete, we consider the example of 
the brick- wall lattice [10] . However, all protocols that we 
propose can be extended to the cases of staggered flux 
lattice [5] and hexagonal lattice [HI E] . 

The scheme for measuring the Zak phase is the follow- 
ing (see Fig. [T]): Initially, a spin-up state with a given 
quasi-momentum ko is prepared. The proposed protocol 
consists of three steps: (1) A 7r/2 pulse is used to create a 
coherent superposition of two spin states, (|t) + /v^; 
(2) Next, opposite forces ±F on the spin-up and spin- 
down are applied, e.g., though a magnetic field gradient. 
It is required that the force F is parallel to some recipro- 
cal lattice vector Gi; (3) After half a period of the Bloch 
oscillations, when the two spins meet in quasi-momentum 
space, another tt/2 pulse is applied and the z-component 
of the spin is measured. During such an evolution, the 
up and down states pick up a geometric contribution, 
equal to the Zak phase, which can be determined from 
the Ramsey phase. 

Bloch oscillations and Zak phase. We start our 
analysis by relating the Ramsey phase measured in the 
setup described above to the Zak phase. Consider a 2D 
lattice with the lattice vectors ai,a2, and the unit cell 
r £ xiBii + X2a.2, where Xi G [0, 1). We denote the two 
primitive reciprocal lattice vectors by Gi, G2. According 
to the Bloch theorem, the eigenfunctions in the nth band 



primitive cell 




FIG. 2: Primitive momentum-space cell of a generic 2D lat- 
tice. The Chern number of the band is related to the winding 
number of the Zak phase across the primitive cell, see Eq. 
(111. The Zak phase can be measured using a combination of 
Bloch oscillations and Ramsey interferometry as described in 
the text. 



can be represented in the following form 

^kn(r) = e'^^Uknir), (1) 

where Ukn is the cell-periodic Bloch function, satisfying 
Wkn(r + a-i) = Wkn(r), i = 1,2. Notice that the Ukn is not 
periodic in the momentum space, but obeys the following 
condition: 



Wk+Giri(,I"j 



(2) 



which originates from the periodicity of the full Bloch 
function V'kn- 

We will assume that during the Bloch oscillations the 
particle motion remains adiabatic, such that the proba- 
bility for the particle to get excited to another band is 
negligible. Following the first 7r/2 pulse, the particles 
are in the state V'kon(r) ® ^^^^^2^- '^^'^ evolution under 
the application of the force ±F (opposite signs for oppo- 
site spins) is described by the time-dependent wave func- 
tion *-t-(r,t) (8) -|J= -I- *;(r, t) (8) The wave functions 
\I'-f(^)(r,t) obey the Schroedinger equation 



(3) 



with the Hamiltonian 



re 

2m 



V{t) being the lattice potential, and Ez the Zeeman en- 
ergy. 

As long as the adiabaticity condition is fulfilled, the 
particles remain within one band, and the solution of 
equation ([3| has the following form: 



vl'^(;)(r,t) = e^«w)(t)^^^(^^,^(r), 



(5) 
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where k±(t) = ko ± ft, f = F/?i, and the phase Ct(^)(*) 
is given by: 

/•k±(t) X ft 

S.ta){t)=i (uk'n|Vk'Uk'n)c?k'-- / e„(k±(t'))di'T 

Jko " Jo 

(6) 

The first term in the above equation describes the ge- 
ometrical phase, while the second and third correspond 
to the dynamical phase, which depends on the speed of 
motion through the band. 

The Ramsey interferometry, performed after half a pe- 
riod of the Bloch oscillations (period is given by T = 
G/|f|), measures the phase difference picked up by the 
two spin species — Using formula (j6]), 

we obtain the Ramsey phase, 

(7) 

where the Zak phase is given by [27]: 

|.ko+G/2 

<PZak = W (wk'ri|Vk'Mk'«)dk' (8) 

"'ko-G/2 

and the dynamical phase and Zeeman phases are given 

by 



(a) 



Ez 



T/2 



T/2 



sign{t')en{ko+it')dt' , (^zc 



^ EzT 

h ■ 

(9) 

For the case of a band structure with symmetric disper- 
sion relation, en(ko + ft') — £„(ko — ft'), the dynamical 
phase vanishes, and the Ramsey interferometry directly 
gives the Zak phase. This is the case for special choices 
of ko and Gi in the experimentally relevant case of the 
brick-wall lattice which we will discuss below. 

Measuring Berry curvature and Chern number 
of a generic band. Let us now turn to the discussion 
of how Ramsey interferometry can be used to determine 
the Berry curvature and the Chern number (and there- 
fore the topological class) of a gapped band; no special 
symmetries are assumed, except for the symmetry of dis- 
persion which guarantees the cancellation of the dynam- 
ical phase, and allows the separation of the Zak phase. 

We choose the primitive cell in quasi-momentum space 
to be a torus defined by k = Kq -I- aiGi -f a2G2, where 
Ui G [0; 1) and Kq is an arbitrary quasi-momentum (as 
shown in Fig. [2]). We notice that the Chern number can- 
not be determined by measuring the Zak phases along 
the four sides of the torus, essentially, because the Zak 
phase is only defined modulo 27r. However, as we now 
discuss, the Chern number C can be related to the wind- 
ing number of the Zak phase across the BZ (see Ref. |2] 
for a closely related discussion in the context of adiabatic 
pumping) . 

We consider an experiment in which the Zak phase is 
measured for torus cycles defined by Gi as a function of 
a2, see Fig. [2] Experimentally, this would be achieved 
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FIG. 3: a) Brick-wall lattice. A and B sites are marked by 
blue and red circles, nearest-neighbor hopping is assumed, b) 
The Brillouin zone of the brick- wall lattice model (blue dashed 
square). The band structure exhibits two Dirac points marked 
by orange circles. Owing to the symmetry of the dispersion, 
it is convenient to measure the Zak phase with initial quasi- 
momentum ko = (fcojO) lying on the x axis, and applying 
a force in the y direction. Measuring the variation of the 
Zak phase as a function of ko, it is possible to (i) measure 
the TT Berry's phase of Dirac particles, (ii) measure the Chern 
number of the bands when they are separated by energy gaps. 



by preparing the initial state ko = Kq + Gi/2 -I- a2G2 
for different values of a2. 

Let us show that the small change of Zak phase as is 
increased by 5a2 is equal to the integral of the Berry cur- 
vature over the rectangle 5S defined by the corresponding 
trajectories (see Fig. [2|. Equivalently, the difference of 
the Zak phases 7 = V3zak(a2 + <^Q!2)-</3zak(a2) is given by 
the Berry's phase that corresponds to the contour 1234. 
It is easiest to see this by choosing a smooth gauge for 
the periodic Bloch function in 5S (this can be done since 
region 5S is small; in general, no smooth gauge can be 
chosen in the whole BZ). The Berry's phase 7 can be 
represented as the sum of the Berry's phases for the four 
sides of the rectangle, 7 — Y^l^ili- Since the sides 2 
and 4 are equivalent (they differ by Gi), but are tra- 
versed in the opposite direction, their contribution van- 
ishes, 72 + 74 = 0. Because we chose the periodic gauge, 
73 + 71 is equal to the difference of the Zak phases for 
trajectories 3 and 1. Thus, the change of the Zak phase is 
related to the Berry's phase, which can be written as an 
integral of the Berry's curvature r2i2, 7 = /^^ (i^fci7i2(k). 
This relation can be conveniently written in terms of a 



uniquely defined quantity z{a2) 



d^kn,2(k) = -iz*{a2)dc,,z{a2)6a2. (10) 



ss 



Summing relation ( 10 ) over different regions, and 



using the definition of the Chern number C = 
Stt /bz '^^'''^^^(k), we then obtain c via the winding of 
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the Zak phase, 



C 



the dispersion (14)). Similarly, it is possible to measure 



2^ 



da2Z*{a2)da2z{a2)- 



(11) 



This relation implies that the interferometric measure- 
ments of the Zak phase across the BZ allow the extrac- 
tion of the Chern number, and can be used to detect the 
topological nature of the bands and topological phase 
transitions. 

Brick-wall lattice: measuring tt Berry's phase 
and detecting topological phase transitions. Our 

approach is also suitable for measuring the tt Berry's 
phase of the massless Dirac fermions in non-trivial lat- 
tices [SHTU]. For definiteness, we consider the brick- wall 
lattice [TD] , illustrated in Fig. |3] 

The particles in the brick-wall lattice with nearest- 
neighbor hoppings t are described by the Hamiltonian 
H = —tJ2{ij) (^l^j (the analysis below can be straight- 
forwardly generalized to the case of anisotropic hopping 
in the x and y direction). Assume that an yl-site coin- 
cides with the origin (this choice is important, because it 
fixes the periodicity condition ([2|, see Fig.|3|. The cell- 
periodic wave function, given by a two-component vec- 
tor {vAk.,VB'k) (corresponding to the amplitudes on the 
sub-lattices A, B), can be found from the Schroedinger 
equation with an effective Hamiltonian 





tt 



t]^ ^ t{2 COS kord + e *"«")= 



(12) 



(13) 



where we have chosen the vector connecting A-site to its 
neighbors to be equal to 1, and the energy of the Bloch 
states is given by e = ±£k, 



£k = t \/4cos2 kj.d + 4 cos fc^-dcos kyd + 1. 



(14) 



The eigenfunctions for the valence and conduction bands 
are given by (uAk,VBk) = ;^(1, ie"'''''). In principle, 
the wave functions are defined up to a gauge transforma- 
tion; here the gauge is fixed such that the wave functions 
satisfy the periodicity condition 

The BZ, illustrated in Fig. [3j is defined by primitive 
reciprocal lattice vectors Gi = {tt,tt) and G2 — (— tt, tt). 
However, measuring the Zak phases for the cycles de- 
fined by G12 is complicated by the lack of corrsponding 
symmetry of the energy dispersion; thus, it would not be 
feasible to cancel the dynamical phase. 

To overcome this difficulty, we propose to measure Zak 
phases for kg lying on the x axis, with the force applied 
in the y direction. This allows for a measurement of the 
Zak phases corresponding to the reciprocal lattice vector 
Hi = Gi -I- G2 (see Fig. [s]), while cancelling the dynam- 
ical phase (as guaranteed by the ky — —ky symmetry of 



the Zak phases for trajectories corresponding to recipro- 
cal vector H2 = Gi — G2, when the initial momentum 
kg lies on the y axis. In this case, the dynamical phase 
cancels out due to kx — >■ —kx symmetry of the dispersion. 

Measuring the Zak phases across the BZ in this case 
allows direct detection of the tt Berry's phase of the Dirac 
fermions. As the initial momentum ko passes either of 
the two Dirac points, situated at (±27r/3, 0), the Zak 
phase must jump by n. 

Furthermore, the method described above can be used 
to determine the topological nature of gapped bands in 
the brick-wall lattice. The winding number of the Zak 
phase measured as a function of initial momentum kg 
on the X-axis gives twice the Chern number (repeating 
the argument above, one can show that Berry's phase 
winding is equal to the integral of Berry's curvature over 
a region — tt < kx-ky < vr, which contains two Bril- 
louin zones). We note that both topologically trivial 
and non-trivial gaps can be introduced by adding ex- 
tra terms in the Hamiltonian (12). The topologically 



trivial gaps arise, e.g., when a staggered potential be- 
tween A and B sites is turned on. In this case, the 
Chern number is zero. The Haldane model, in which 
next-nearest-neighbor hoppings with non-trivial phases 
±77 are turned on |35) . provides a realization of a non- 
trivial band, with the Chern number ±1. Importantly, 
both models still preserve the symmetry of dispersion 
with respect to kx — >■ —kx,ky —ky (the Haldane model 
obeys this condition when 77 = ±7r/2). Thus, the prob- 
lem of dynamical phases does not occur when the Zak 
phases (and Chern numbers) are measured. 

Protocols insensitive to fluctuating magnetic 
fields. Now we would like to address an important exper- 
imental concern related to the presence of slowly fluctu- 
ating magnetic fields. The fluctuations of the magnetic 
field induce shot-to-shot variations in Ez and ifzeemam 
complicating the reliable extraction of the Zak phase. 
We propose various protocols which are insensitive to the 
magnetic noise. Here we briefly mention the main ideas, 
and provide a detailed description in Ref. [36j. 

The flrst protocol allows measurements of the varia- 
tions of the Zak phase across the BZ (and in fact, only the 
difference of the Zak phase is physically meaningful, while 
the Zak phase itself depends on the choice of the real- 
space unit cell, see, e.g. [31]). The idea is to prepare the 
system in two or more different quasi-momentum states, 
and to carry out the sequence described above for each 
initial state. Assuming that the dispersion is symmet- 
ric, the difference of the phases extracted from Ramsey 
fringes will be equal to the difference of the Zak phases. 
This protocol therefore allows one to obtain the Berry 
curvature and the Chern number of the band. It is in- 
sensitive to the shot-to-shot fluctuations of the magnetic 
field because the Zeeman phases picked up by different 
quasi-momentum states are equal. A more detailed dis- 
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cussion of this protocol for the brick-wah lattice can be 
found in Ref . [36j . Another approach is to design proto- 
cols that combine Bloch oscillations with spin echo se- 
quence [36] . Such protocols are naturally insensitive to 
the fluctuating magnetic fields. 

Summary. In conclusion, we have presented an ap- 
proach to studying topological properties of 2D optical 
lattices. This approach allows one to measure the Berry's 
phase of the Dirac fermions, as well as the local Berry 
curvature and the Chern number of the band. The latter 
measurement can be used to study the topological struc- 
ture of the bands in optical lattices and to detect the 
topological phase transitions that occur as a function of 
the lattice parameters. 
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SUPPLEMENTARY ONLINE MATERIAL FOR "INTERFEROMETRIC APPROACH TO MEASURING 

BAND TOPOLOGY IN 2D OPTICAL LATTICES" 

Here, we provide a detailed description of the protocols which allow one to eliminate the effect of fluctuating 
magnetic fields. We illustrate the ideas by using the experimentally relevant example of brick-wall lattice, noting that 
extensions to the case of flux lattice 8 and honeycomb lattice are straightforward. 

The first protocol (which was discussed in the main text in some detail) is illustrated in Fig. |4] One can envision that 
this protocol can be implemented using either bosons or fermions. For the case of bosons, one would prepare two clouds 
of atoms with different quasi-momenta kx = k-mim kx = fcmax and ky = 0. We cannot predict the actual Ramsey/Bloch 
phases in a single run of the experiment due to magnetic field fluctuations which lead to fluctuations of the Zeeman 
contribution. However, the difference of the phases should be n when the two quasi-momentum states lie on the 
opposite sides of the Dirac point, as in Fig. [4] Note that particle crossing the Dirac point do not have a well defined 
Ramsey/Bloch phase due to the breakdown of adiabaticity condition. Alternatively, an ensemble of cold fermionic 
atoms occupying a finite region in quasi-momentum space, is loaded into a brick-wall lattice. The Ramsey/Bloch 
phases can be measured for each momentum k separately. Thus we can compare accumulated Ramsey/Bloch phases 
for different initial momenta. This allows measurement of the tt Berry's phase of Dirac fermions, and, more generally, 
the Berry's curvature. 

Another approach is to combine Bloch oscillations with spin-echo sequence. Such protocols are naturally insensitive 
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FIG. 4: A protocol in which a variation of Zak phase across BZ in measured, for the case of a brick-wall lattice, (a) An 
ensemble of atoms with different quasi-momenta is loaded into an optical lattice (grey cloud). The Ramsey/Bloch sequence 
allows measurement of the Zak phase for quasi-momenta in the interval fcmin < kx < fcmax. (b) The Zak phase exhibits a tt 
jump when crosses the Dirac point (situated at kx = 2n/3), which provides a signature of the n Berry's phase of Dirac 
fermions. 




FIG. 5: (a) Spin-echo protocol for measuring the Berry's phase of the Dirac fermions and Chern number of the band, (b) 
Modified spin-echo protocol for measuring the Berry's phase of the Dirac fermions in brick-wall lattice. 



to fluctuating magnetic flelds, because each atom spends equal time being in the spin-up and spin-down state. 
Therefore, the Zeeman phases picked up by all atoms are equal and cancel each other out. 

One possible spin-echo-type protocol, illustrated in Fig. [sja), consists of three steps: (1) A n/2 pulse, combined 
with half a period of Bloch oscillations, similar to the sequence for measuring the Zak phase discussed in the main 
text. (2) Then, both spin components are quickly moved in the y direction (applied forces are equal). This is followed 
by the tt pulse, which flips spins up and down. (3) Another half a period of Bloch oscillations is performed. (Notice 
that the direction of the forces apphed to pseudo-spin up and down is the same as in step (1).) After that, a 7r/2 
pulse is applied and the accumulated phase is measured. When the trajectory of the particle encloses the Dirac point 
(as shown in Fig. [sjja)), the accumulated phase is equal to tt (the Berry's phase of Dirac fermions). We emphasize 
that this protocol is not limited to the case of brick-wall lattice, and in general allows the measurement of the band's 
Berry's curvature and the Chern number. 

Another possible protocol (which we call modified spin-echo) uses the experimental ability to tune the band structure 
in situ [TU]. The idea of the protocol, illustrated in Fig. [sjb), is not to move atoms in the same direction (step (2) 
described above), but rather to move Dirac points by dynamically changing the band structure. The protocol consists 
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of three steps: (1) A 7r/2 pulse followed by half a period of Bloch oscillations. (2) A tt pulse. (3) Half a period of Bloch 
oscillations, followed by 7r/2 pulse and measurement of the Ramsey phase. During steps (l),(2),and (3) parameters 
of the optical lattice are changed adiabatically such that atoms return on the other side of the Dirac point. This 
measurement gives the tt Berry's phase of the Dirac fermions. 



